home *** CD-ROM | disk | FTP | other *** search
/ Collection of Tools & Utilities / Collection of Tools and Utilities.iso / fortran / pcrng.zip / RMUNI.ASM < prev    next >
Assembly Source File  |  1990-06-29  |  15KB  |  367 lines
  1.               .8087
  2.               page ,132
  3. ;
  4. ;***********************************************************************
  5. ;* Here are data declarations for use by RINIT, IVNI, IUNI, VNI, UNI.  *
  6. ;* All routines are written for use with IBM Fortran/2 Compiler.       *
  7. ;*                                                                     *
  8. ;*   Authors: G. Marsaglia, B. Narasimhan and Arif Zaman               *
  9. ;*                 Supercomputer Computations Research Institute       *
  10. ;*                         and                                         *
  11. ;*                 Department of Statistics                            *
  12. ;*                 Florida State University                            *
  13. ;*                 Tallahassee, Fl 32306-3033.                         *
  14. ;***********************************************************************
  15. ;
  16.               .MODEL     small                        ;Use small memory model
  17.               TITLE      'A Random Number Generator for PC's.'
  18. ;
  19. ;---------------Equates-----------------
  20. ;
  21. arg_offset    equ            6                        ;Offset of argument.
  22. c_low         equ       0cbb1h
  23. c_high        equ          74h                        ;C = 7654321 initially.
  24. cons_low      equ        55e5h
  25. cons_high     equ        159ah                        ;Constant = 362436069.
  26. no_of_bits    equ           32                        
  27. top_of_list   equ          168
  28. init_i_val    equ  top_of_list
  29. init_j_val    equ           84
  30. ;
  31. datseg        segment     para 'F@DATA'
  32. ;
  33. ; Table for Generator (Default values).
  34. ;
  35. unitabl       dd         0F83CEE7Bh
  36.               dd         0A83E5AD3h
  37.               dd          036200BBh
  38.               dd         0FA5764F6h
  39.               dd         0A13CBFC4h
  40.               dd          565A191Eh
  41.               dd          14D4CCFBh
  42.               dd          7F5AD22Ch
  43.               dd          03528F2Eh
  44.               dd         0E81E32DDh
  45.               dd          71C47276h
  46.               dd         0AA0F8045h
  47.               dd          3C3F1C78h
  48.               dd         0E8CE101Dh
  49.               dd         0CCC12691h
  50.               dd          47196DBFh
  51.               dd          074A6DFFh
  52.               dd          03FB675Eh
  53.               dd          60436236h
  54.               dd          072F0247h
  55.               dd          0C7E9185h
  56.               dd         0CE579EAEh
  57.               dd          01864E96h
  58.               dd         0A526C5C6h
  59.               dd         0F582EB0Dh
  60.               dd         0AFE827F4h
  61.               dd          55CED836h
  62.               dd          26124C49h
  63.               dd          7049AEE1h
  64.               dd          49552795h
  65.               dd         0D1602ED6h
  66.               dd          051C65CEh
  67.               dd         0AEF3CC37h
  68.               dd         0A83880ABh
  69.               dd         0AE7EE06Ah
  70.               dd         0D64D988Bh
  71.               dd          0A96BC73h
  72.               dd         0ECEF4297h
  73.               dd          6C18300Dh
  74.               dd          22F3A897h
  75.               dd         0B4D760ADh
  76.               dd         0AC838383h
  77.               dd         0FD04E68Fh
  78. ;
  79. i             dw               164
  80. j             dw                80                    ;Indices I and J.
  81. c             dd         0EADA75CCh                   ;Constant C.
  82. mod169        db               169                  
  83. mod179        db               179                  
  84. multiplier    db                53
  85. carry         db                 1                    ;Carry bit (Default).
  86. stored_val    dd          122a70c3h
  87. datseg        ends
  88. ;
  89. ;***********************************************************************
  90. ;* Subroutine RINIT takes an 8-digit number as argument. Use as        *
  91. ;*                  CALL RINIT(I)                                      *
  92. ;*                                                                     *
  93. ;*  Purpose: Seeds the table and prepares pointers for use by other    *
  94. ;*           routines.                                                 *
  95. ;*                                                                     *
  96. ;* This is set up for use with IBM Fortran/2 Compiler.                 *
  97. ;* -----Note that a 8087 co-processor is assumed to be present.------- *
  98. ;*                                                                     *
  99. ;*   Authors: G. Marsaglia, B. Narasimhan and A. Zaman                 *
  100. ;*            Supercomputer Computations Research Institute            *
  101. ;*                    and                                              *
  102. ;*            Department of Statistics                                 *
  103. ;*            Florida State University                                 *
  104. ;*            Tallahassee, Fl 32306-3033.                              *
  105. ;***********************************************************************
  106. ;
  107. rseg          segment    'CODE'
  108. ;
  109.               assume cs:rseg, ds:datseg
  110. rinit         proc       Far
  111.               public     rinit
  112. ;
  113. ; Usual IBM Fortran/2 stuff....
  114. ;
  115.               push       bp                           ;Save caller's BP.
  116.               mov        bp,sp                        ;Ready to address.
  117. ;
  118. ; No traceback needed.
  119. ;
  120.               mov        ax,datseg                    ;DATA segment.
  121.               mov        ds,ax
  122. ;
  123. ; Load parameter and strip two digits at a time to get four seed values.
  124. ;
  125.               les        si,DWORD PTR[bp].arg_offset
  126.               mov        ax,Word Ptr es:[si]
  127.               mov        dx,Word Ptr es:[si]+2
  128.               mov        cx,10000
  129.               div        cx 
  130.               mov        cl,100
  131.               div        cl
  132.               mov        bx,ax                        ;BH = K, BL = L.
  133.               mov        ax,dx
  134.               div        cl                    
  135.               mov        dl,ah                
  136.               mov        dh,al                        ;DH = J, DL = I.        
  137. ;
  138. ; Add one to all to make sure they are non-zero.
  139. ;
  140.               inc        dl
  141.               inc        dh
  142.               inc        bh
  143.               inc        bl
  144. ;
  145. ; DL = I, DH = J, BL = L, BH = K.
  146. ;
  147.               xor        si,si                        ;Table Index.
  148. labl01:
  149. ;
  150. ; Result will be in [BP,DI].
  151. ;
  152.               xor        di,di                        
  153.               xor        bp,bp
  154.               mov        cx,no_of_bits                ;Bit counter.
  155. labl02:
  156.               shl        di,1              
  157.               rcl        bp,1                         ;Result *= 2.
  158.               mov        al,dh                        ;AL = J.
  159.               mul        bh                           ;AX = J * K.
  160.               div        mod179
  161.               mov        al,ah                        ;AL = AX mod 179.
  162.               mul        dl                           ;AX = AX * I.
  163.               div        mod179
  164.               mov        dl,dh                        ;I = J.
  165.               mov        dh,bh                        ;J = K.
  166.               mov        bh,ah                        ;K = new K.
  167. ;
  168.               mov        al,bl                        ;AL = L.
  169.               mul        multiplier
  170.               inc        ax
  171.               div        mod169
  172.               mov        bl,ah                        ;L = 53*L+1 mod 169.
  173. ;
  174.               mov        al,bl                        ;AL = new L.
  175.               mul        bh                           ;AX = L * K.
  176.               and        ax,63                        ;Mod 64.
  177.               cmp        ax,32
  178.               jl         labl03
  179.               inc        di                           ;Set low bit.
  180. labl03:
  181.               loop       labl02          
  182. ;
  183.               mov        WORD PTR unitabl[si],di       
  184.               mov        WORD PTR unitabl[si+2],bp    ;Save value.
  185.               add        si,4                         ;Bump pointer.
  186.               cmp        si,172                       ;Seeded 86 entries?
  187.               jne        labl01
  188. ;
  189. ; Set indices.
  190. ;
  191.               mov        i,init_i_val
  192.               mov        j,init_j_val
  193. ;
  194. ; Set ca